Metabolomic changes in crown of alfalfa (Medicago sativa L.) during de-acclimation

Alfalfa is a high-quality forage legume species that is widely cultivated at high latitudes worldwide. However, a decrease in cold tolerance in early spring seriously affects regrowth and persistence of alfalfa. There has been limited research on the metabolomic changes that occur during de-acclimation. In this study, a liquid chromatography–mass spectrometry system was used to compare the metabolites in two alfalfa cultivars during a simulated overwintering treatment. In four pairwise comparisons, 367 differential metabolites were identified, of which 31 were annotated according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Many of these metabolites were peptides, carbohydrates, and lipids. At the subclass level, 17 major pathways were revealed to be significantly enriched (P < 0.05). The main differential metabolites included amino acids, peptides and analogs, carbohydrates, and glycerol phosphocholines. A metabolomic analysis showed that the up-regulation of unsaturated fatty acids and amino acids as well as the enhancement of the related metabolic pathways might be an effective strategy for increasing alfalfa cold tolerance. Furthermore, glycerophospholipid metabolism affects alfalfa cold tolerance in early spring. Study results provide new insights about the changes in alfalfa metabolites that occur during de-acclimation, with potential implications for the selection and breeding of cold-tolerant cultivars.

In temperate climates, seasonal changes of alfalfa (Medicago sativa L.) include cold acclimation in autumn, dormancy in winter, and de-acclimation and regrowth in spring. For optimal survival, alfalfa acquires an adequate level of freezing tolerance, which is maintained for several months, and then de-acclimates after early spring freeze-thaw cycles. In northern China, because of the extended periods of sub-zero air and soil temperatures, only cold tolerant cultivars can be grown. These cultivars can tolerate extremely low temperatures during winter, but they can also adapt to temperature fluctuations in spring 1,2 .
A series of physiological, biochemical, and molecular changes induced by photoperiod and temperature maximize the freezing tolerance of plants during acclimation [3][4][5] . Earlier research indicated that the hardiness of perennial ryegrass (Lolium perenne L.) 6 and asparagus (Asparagus officinalis L.) 7,8 is negatively associated with increasing water contents. A decrease in the water content can lower the freezing point of cells, whereas an increase in the water content may result in mechanical damages due to extracellular freezing as well as an increase in the ice propagation rate in tissues [7][8][9] . Furthermore, the "energy-competitive hypothesis" states that the allocation of substances during acclimation will affect the freezing tolerance of plants, with a decrease in metabolizable energy consumption and an increase in energy production conducive to low-temperature stress tolerance when the water content is low 5 . The soluble sugar content in alfalfa root crowns reportedly increases several-fold in a few hours under low-temperature conditions 10 . Additionally, the superior winter hardiness of some cultivars is related to the early initiation of the accumulation of soluble sugars in autumn and high sugar levels 11,12 . Previous studies revealed that the differences between non-hardy and winter-hardy alfalfa cultivars are related to the capacity of plants to accumulate starch, sucrose, glucose, stachyose, raffinose, maltose, and fructose 8,[10][11][12][13] . These soluble sugars prevent plants from being damaged by an exposure to low temperatures by maintaining the cell membrane integrity 8 , transmitting cold signals 14 , and regulating the expression of lightinduced genes 3 during acclimation. The plasma membrane is the primary site of freezing damage, the extent of which is closely associated with plasma membrane characteristics 15,16 . Thus, the plasma membrane lipid composition is an important factor affecting freezing tolerance. The molecular species composition of plasma membrane lipids has been investigated in Arabidopsis (Arabidopsis thaliana L.) 15 , winter wheat (Triticum aestivum L.) 17 ,

Materials and methods
Plant materials. This study was conducted at the Institute of Animal Science, Chinese Academy of Agricultural Sciences (Beijing, China) from August 2020 to February 2021. The two alfalfa cultivars selected for this study were 'WL440' (fall dormancy score: approximately 6) and 'ZhaoDong' (ZD; fall dormancy score: approximately 2), which are commonly grown in northern China. These fall dormancy scores are representative of the alfalfa cultivars grown in the region. The WL440 seeds were provided by Beijing Zhengdao Seed Industry Co., Ltd. (Beijing, China), whereas the ZD seeds were obtained from Gansu Agricultural University (Lanzhou, China). Seeds were disinfected in sodium hypochlorite (1% NaClO) for 30 min and then washed five times with deionized water. Seeds similar in size were selected and germinated for 48 h in a Petri dish (diameter of 9 cm) at 25 °C. Eight germinated seeds were transferred into a polyvinyl chloride (PVC) pipe (inner diameter of 11 cm and length of 25 cm). Each PVC pipe was filled with 600 g sterilized dry mixture of sandy soil and nutrient soil (4:1, v/v). The dry mixture consisted of peat moss and lime, with 1.6% total nitrogen, 0.1% P 2 O 5 , 0.2% K 2 O (N:P:K = 14:10:18), 91% organic matter, conductivity of 0.9 dS m −1 . The water-holding capacity of the mixture was 38.35%. Four plants that were about 20 cm tall were retained in each PVC pipe 2 weeks later. They were then allowed to grow for an additional 5 months, during which mowing was performed on September 30, November 10, and December 20, 2020, with 10 cm of stubble left after each mowing. The pipes were weighed every 4 days to check that the soil moisture content was maintained between 60 and 65% of water-holding capacity. Weeds and pests were removed regularly. The use of plants in the present study complies with international, national and/or institutional guidelines.
Experimental design and sampling. The study was completed using a randomized complete block design involving two alfalfa cultivars and four experimental phases ( Fig. 1) (20 °C); 350 µmol·m −2 ·s −1 (photosynthetic photon flux density); 65% relative humidity. Phase 1 sampling occurred on January 20, 2021. The remaining materials were transferred to the LRH-200-GD low-temperature illumination incubator (Taihong Medical Instruments, Guangdong, China) for the subsequent analyses. The initial incubation conditions were as follows: 14 h day (25 °C)/10 h night (20 °C); 350 µmol·m −2 ·s −1 (photosynthetic photon flux density). The temperature was decreased to 5 °C (day)/0 °C (night) at a rate of 2 °C·d −1 and the light intensity was decreased to 150 µmol·m −2 ·s −1 at a rate of 20 µmol·m −2 ·s −1 ·d −1 to simulate the cold acclimation environment for alfalfa in autumn. Phase 2 sampling was performed after an additional 72 h of cold acclimation. The temperature was decreased further to − 5 °C (day)/ − 10 °C (night). Phase 3 sampling was conducted after a 72 h incubation. The temperature was then increased to 10 °C (day)/5 °C (night) at a rate of 2 °C·d −1 . Phase 4 sampling was performed on February 21, 2021. The aboveground and belowground plant parts were separated and the root was carefully rinsed by hand with distilled water. The root crown (about 5 cm belowground) was divided into two parts. One part was used for measuring electrical conductivity, whereas the other part was stored at − 80 °C for the metabolomic analysis. The remainder of the root was measured to calculate the biomass. The aboveground and belowground parts were weighed after being placed in an oven at 65 °C for 48 h. Dry weights were recorded as the aboveground and belowground biomass.
Determination of the semi-lethal temperature. The temperature at which the relative permeability of intracellular ions reaches 50% (LT 50 ) was used to represent the cold hardiness of alfalfa 38 . The root-crown portion of plants were sliced into nine pieces (2-3 mm each) and placed in separate 2 mL centrifuge tubes.
Tubes were kept at 4 °C prior to testing. The subsequent freezing test was conducted using the ZX-5C constanttemperature circulator (Zhixin Instrument, Shanghai, China). Sample tubes were incubated in alcohol for 1. Electrical conductivity (i.e., EL 1 ) was measured using a FE38 conductivity meter (Mettler, Shanghai, China). Samples were autoclaved at 121 °C for 30 min and electrical conductivity was remeasured (i.e., EL 2 ). Electrical conductivity of the deionized water was also measured (i.e., EL). Relative electrolyte leakage was calculated according to Eq. (1) and the semi-lethal temperature was calculated using logistic Eq. (2), in which x is the freezing temperature, y is the relative electrolyte leakage, and A, B, and k are constants:

Metabolite extraction.
To determine the changes in metabolite contents during de-acclimation, phase 3 and 4 samples (i.e., de-acclimation stages) were analyzed. Phase 3 and 4 samples of ZD and WL440 were designated ZD_1, ZD_2 and WL440_1, WL440_2. Frozen crown samples were ground to a fine powder in liquid nitrogen. 50 mg was transferred to a 2 mL centrifuge tube and a grinding ball (diameter of 6 mm) was added. 400 μL extraction solution was added composed of methanol and water (4:1, v/v) and 0.02 mg mL −1 L-2-chlorophenylalanine. Samples were ground for 6 min at − 10 °C using a grinder (50 Hz). An ultrasonic extraction (40 kHz) was performed for 30 min at 5 °C and samples were transferred to a centrifuge tube and incubated at − 20 °C for 30 min. The mixture was centrifuged at 13,000 g at 4 °C for 15 min. 200 μL supernatant was then added to an LC-MS sample vial. 20 μL supernatant aliquots from each sample were mixed to prepare quality control samples. To monitor the stability and repeatability of the instrumental analysis, the quality control samples (inserted after every eight samples) were analyzed along with an equal volume of each sample.

Liquid chromatography-mass spectrometry (LC-MS) analyses.
The LC-MS analyses were performed using ultra-high performance liquid chromatography (UHPLC)-Q Extactive HF-X system (Thermo Fisher Scientific, Waltham, MA, USA). Chromatographic separation was completed using the ACQUITY UPLC HSS T3 column (100 mm × 2.1 mm i.d., 1.8 μm; Waters, Milford, USA). Mobile phase A was 95% water and 5% acetonitrile (containing 0.1% formic acid). Mobile phase B was 47.5% acetonitrile, 47.5% isopropanol, and 5% water (containing 0.1% formic acid). The injection volume was 2 μL. The column temperature was 40 °C. Details regarding the mobile phase elution gradient are provided in Table S1. Samples were ionized (via electrospray) and mass spectrometry signals were collected in positive and negative ion scanning modes. The specific parameters are provided in Table S1.
Bioinformatics and data analyses. Data were analyzed using the metabolome data processing software Progenesis QI (Waters Corporation) (i.e., for baseline filtration, peak identification, integration, retention time correction, and peak alignment). The data matrix containing retention time, mass-to-charge ratio, and peak intensity were identified. MS and MS/MS data were used to screen for matches in the metabolite databases. The The data matrix was converted using SIMCA-P software (version 14; Umetrics, Umea, Sweden) to perform a principal component analysis (PCA), a partial least squares discriminant analysis (PLS-DA), and an orthogonal partial least squares discriminant analysis (OPLS-DA). The criteria used for identifying significant differential metabolites were as follows: VIP (Variable Importance in Projection) > 1 and P < 0.05. Differential metabolites were analyzed qualitatively and the compound identification numbers (CIDs) were acquired using the Human Metabolome Database (https:// hmdb. ca/). Differential metabolites were further annotated using the KEGG compound database (http:// www. kegg. jp/ kegg/ compo und/) [39][40][41] . Metabolites were annotated according to the CIDs and were mapped using the KEGG pathway database (https:// www. genome. jp/ kegg/ pathw ay. html). Metabolites were classified, prominent metabolic pathways were identified, and metabolite enrichment was analyzed.
Analysis of variance, using SPSS 20.0 (SPSS Inc., Chicago, IL, USA), was conducted for biomass and LT 50 data for the cultivars and phases. Multiple range tests were performed using least significant differences. The threshold for determining the significance of any differences was P < 0.05.

Results
Biomass and semi-lethal temperature. For both cultivars, the aboveground biomass increased from phase 1 to 4, but no significant effect was detected in the within-cultivar comparisons (Fig. 2a,b). The growth rate of WL440 was greater than that of ZD in phase 4. The aboveground biomass of ZD was significantly greater than that of WL440 in phases 1 and 2. There was no significant difference in the belowground biomass of the two cultivars in any phase. From phase 1 to 4, the belowground biomass first increased and then decreased. The belowground biomass was highest in phase 2 for ZD and in phase 3 for WL440. The crown water content decreased from phase 1 to 3 and then was stable until phase 4. The crown water content did not differ significantly between the two cultivars, but there were significant differences among the four phases (Fig. 2c). There was no significant difference in the semi-lethal temperature between the two cultivars in phase 1, but there were significant differences from phase 2 to 4 (Fig. 2d). The difference in the semi-lethal temperature between the two cultivars was greatest in phase 2, after which it decreased in phase 3 and then increased in phase 4.
Classification of differential metabolites. PCA data (Fig. 4) show that the differential metabolites were well separated between phases and cultivars and the variation within groups was small. Using VIP > 1.0 and P < 0.05 as criteria, 367 differential metabolites were identified (Table S3) and classified as follows: lipids and lipid-like molecules (133), organic acids and derivatives (75), organic oxygen compounds (44), phenylpropanoids and polyketides (42), organic heterocyclic compounds (39), benzenoids (16), nucleosides, nucleotides, and analogs (7), lignans, neolignans, and related compounds (6), organic nitrogen compounds (2), alkaloids and derivatives (1), homogeneous non-metal compounds (1), and hydrocarbon (1). Volcano maps of the differential metabolites in different pairwise comparisons are presented in Fig. 5. There were 27 up-regulated and 80 down-regulated differential metabolites in the WL440_2 vs WL440_1 comparison, whereas there were 42 up-regulated and 32 down-regulated differential metabolites in the ZD_2 vs ZD_1 comparison. The within-cultivar comparisons revealed that in WL440, the number of down-regulated metabolites  www.nature.com/scientificreports/ increased when the temperature increased. Additionally, 25 up-regulated and 66 down-regulated differential metabolites were identified in the ZD_1 vs WL440_1 comparison, whereas 54 up-regulated and 27 down-regulated differential metabolites were identified in the ZD_2 vs WL440_2 comparison (Table 1). Specific details regarding the differential metabolites identified in the four comparison groups (e.g., classifications and fold changes) are provided in Table S4. The differing metabolites in the four pairwise comparisons (Fig. 6) were divided into 33 classes and 58 subclasses (details are provided in Table S5). Lipids and lipid-like molecules, organic acids and derivatives, and organic oxygen compounds were the most active metabolites according to the proportions of all classes of differential metabolites. At the subclass level, the differential metabolites were primarily organooxygen compounds, Figure 5. Volcano maps of differential metabolites in four pairwise comparisons: (a) WL440_2 vs WL440_1; (b) ZD_2 vs ZD_1; (c) ZD_1 vs WL440_1; (d) ZD_2 vs WL440_2. Note: Volcano map is used to show the relative content difference of metabolites between the two groups of samples and the statistical significance of the difference. Each point in the volcano map represents a metabolite, where red points represent up-regulated differential metabolites, gray points for non-significantly different metabolites, and green points for downregulated differential metabolites. The abscissa represents the logarithmic value (log 2 FC) of the multiple difference of the relative content of a metabolite in the two groups of samples. The larger the absolute value of the abscissa, the greater the relative content difference of the metabolite between the two groups of samples. The ordinate represents the VIP value, and the larger the ordinate value, the more significant difference is and the more reliable the differential metabolites obtained by screening. www.nature.com/scientificreports/ prenol lipids, fatty acyls, carboxylic acids and derivatives, and glycerophospholipids. Carboxylic acids and derivatives, organooxygen compounds, and prenol lipids were the main differential metabolites that were up-regulated in the WL440_2 vs WL440_1, ZD_2 vs ZD_1, and ZD_2 vs WL440_2 comparisons, whereas carboxylic acids and derivatives, prenol lipids, and glycerophospholipids were the main up-regulated differential metabolites in the ZD_1 vs WL440_1 comparison. Carboxylic acids and derivatives and prenol lipids were identified as downregulated differential metabolites in all pairwise comparisons. The other identified down-regulated differential metabolites included fatty acyls in the WL440_2 vs WL440_1 and ZD_2 vs ZD_1 comparisons, glycerophospholipids in the ZD_2 vs ZD_1 and ZD_2 vs WL440_2 comparisons, and organooxygen compounds in the ZD_1 vs WL440_1 and ZD_2 vs WL440_2 comparisons.  Figure 6. Classification of differential metabolites in four pairwise comparisons. Note: The abscissa represents the names of the 58 subclasses level of differential metabolites in the four pairwise comparisons.  (Fig. 7). These metabolites were divided into seven groups: peptides (14), carbohydrates (5), lipids (4), nucleic acids (3), hormones and transmitters (3), vitamins and cofactors (1), and organic acids (1). According to the KEGG pathway database, 17 major pathways were significantly enriched (P < 0.05) among the differential metabolites (Fig. 8). The pathways that were identified as enriched, but not significantly are listed in Table S6. At the subclass level, the metabolites associated with the 17 major pathways included amino acids, peptides and analogs, carbohydrates and carbohydrate conjugates, and glycerophosphocholines ( Table 2).  www.nature.com/scientificreports/ In conclusion, carbohydrates, amino acids and lipids are the important differential metabolites responsible for the difference in cold resistance between ZD and WL440, and these metabolites are of great significance for breeding new cold-resistant alfalfa cultivars.

Discussion
Although there is no consensus regarding the relationship between the plant tissue water content and cold hardiness, most researchers believe that cold hardiness is negatively associated with increasing water content due to lower water content can decreased osmotic pressure 5,7,8,42 . In this study, an analysis of acclimation revealed that the freezing tolerance of alfalfa was highest in phase 3 because of its decreased water content and freezing conditions (Fig. 2c,d). Additionally, cellular dehydration reportedly decreases metabolic activity, growth, and energy consumption 5 . De-acclimation and regenerative growth are related to tissue rehydration. Compared with tissues with low water contents, tissues with high water contents are more susceptible to mechanical damage caused by extracellular freezing and have a higher ice propagation rate 23 . Table 2. Changes of the differential metabolites in four pairwise comparisons that involved in 17 pathways identified as being significantly enriched in the KEGG analysis. "↓" represents a metabolite was downregulated in two comparison groups, while "↑" represents a upregulation. "-" indicates that the metabolite or a significant difference was not detected.

Terpene glycosides
Geniposidic Acid www.nature.com/scientificreports/ Alfalfa survival during a cold winter can be attributed to several causes 13,22,30 : (1) decreases in the temperature and photoperiod induce acclimation in alfalfa; (2) a continuously decreasing temperature causes alfalfa to almost stop growing, acquire peak cold hardiness, and survive the midwinter; (3) an increase in the temperature in early spring decreases the cold hardiness of alfalfa and induces regeneration. Moreover, temperature fluctuations may induce re-acclimation and de-acclimation 1,22 . Plants have evolved various strategies to minimize cellular damage in response to low-temperature stress. Examples include osmotic regulation, synthesis of soluble proteins, changes in cell membrane components, and hormonal regulation 1 .
In the current study, there was no significant difference in the cold hardiness of ZD and WL440 in phase 1, but a significant difference was detected in phase 2, possibly due to the differential metabolism of soluble sugars and amino acids 13,43 . Soluble sugars are sensitive to temperature changes. An earlier study by Bertrand 10 demonstrated that the soluble sugar content in the root-crown region of alfalfa increases several-fold after only 8 h of a low-temperature treatment. Similarly, several recent studies have confirmed the existence of photoperiodmeasuring mechanisms in the plant metabolic network that coordinates seasonal developmental programs to mitigate damages from less predictable stresses 3,4,44 . Liu 4 reported that in Arabidopsis, changes in the cellular starch and sucrose contents can regulate the expression of photoperiod-related genes, thereby maintaining the diurnal energy level, which is crucial for plants to adapt to the relatively long and cold nights in winter. Amino acids serve as a nitrogen source for plants, but they are also important for adaptations to low-temperature stress 45,46 . Amino acids accumulate in the alfalfa root-crown region at low temperatures and specific amino acids can delay protein degradation and maintain the structural stability of the plasma membrane 47-50 . In phase 3, there was a minor difference between the two cultivars, suggesting that sufficient cold acclimation will decrease the differences in the sensitivity to low temperatures among cultivars. Research regarding climate changes due to global warming has generated evidence that temperature fluctuations increase the risk of overwintering failure compared with the effects of an exposure to consistent low temperatures 51,52 .
Temperature increases result in changes in a series of plant metabolic pathways, which lead to the rapid loss of cold hardiness. In the current study, fucose content increased after de-acclimation. This is contrary to earlier findings 13 . As a signaling molecule and protective agent, fucose regulates plant growth and development and the mechanism mediating responses to environmental stimuli. In addition to protecting proteins and cell membranes at low temperatures, fucose also co-regulates the plant heat stress response mechanism with hormones when temperature increases [53][54][55] . The composition and fluidity of the plasma membrane are key factors influencing plant sensitivity to low temperatures. At low temperatures, the plasma membrane transforms from a relatively fluid structure to a solid gel, in which lipids are tightly packed and highly ordered. The solid gel membrane is highly permeable and is easily ruptured 54 . The decrease in glycerophospholipid metabolic activity during deacclimation is believed to be related to the loss of cold hardiness. Glycerophosphate lipids, which are phospholipids comprising a phosphate and two fatty acids or fatty alcohols, are the main plasma membrane components. The comparison of different phases in this study revealed decreases in glycerophosphate, glycerophosphocholine, lysophosphatidylcholine (16:0), and lysoPC (15:0) content during de-acclimation. The decrease in unsaturated fatty acid content accelerated the loss of cold hardiness. Unsaturated fatty acids are useful for maintaining plasma membrane fluidity and cellular functions under low-temperature stress conditions, likely because of their cisdouble bonds 13,15,54 . Miki 56 conducted a proteomics-based analysis of Arabidopsis, which confirmed that changes in amino acid contents in plants are closely related to temperature. During acclimation and de-acclimation, hundreds of proteins are modified, among which transporters are the most reduced. On the basis of an analysis of the circadian clock and energy, Bonnot 44 reported that in response to temperature limitation, organisms selectively translate mRNA to respond to environmental changes. According to a cluster analysis, most of the changes to proteins during acclimation and de-acclimation were reversible. Additionally, the changes during de-acclimation were faster than those during acclimation. In phase 4, tyrosine, succinic acid, and methionine decreased significantly in WL440, but were not significantly changed in ZD, this may help explain the observed differences between these cultivars. These findings were in accordance with the results of a previous study 57 , in which functional amino acids were observed to protect cells from cold damage.
The metabolic pathways that differed between ZD_2 and WL440_2 were mainly related to purine metabolism and plant hormone signal transduction (Fig. 8). Abscisic acid (ABA) is a sesquiterpenoid hormone that affects plant growth and development, while also regulating plant responses to adverse environmental conditions. Furthermore, ABA-dependent signaling is one of the pathways responsive to low-temperature stress [58][59][60] . An exposure to low temperatures leads to an increase in ABA content. The accumulated ABA binds to receptors, which then interact with PP2C (type 2C protein phosphatase) to inhibit the binding of PP2C to SnRK2 (SNFIrelated protein kinase 2). As a result, SnRK2 transcription factors are phosphorylated, which can activate the expression of ABA-responsive genes and increase cold hardiness 61,62 .
The mechanisms underlying alfalfa plant responses to cold conditions are primarily associated with changes in amino acids, lipid molecules, carbohydrates, and their related metabolic pathways. At different stages, various metabolites contribute to alfalfa cold hardiness. For example, alfalfa adapts to low-temperature environments via an increase in the abundance of amino acids (e.g., aspartic acid), fatty acids (e.g., arachidonic acid and norlinolenic acid), glycerophospholipids (e.g., lysophosphatidylcholine and phosphatidyl ethanolamine), and carbohydrates (e.g., glucose, maltose, and raffinose) in the root-crown region during phase 2 13 . However, in phase 3, the contents of most of these metabolites are no longer increased, and may be decreased (with the exception of raffinose). Cystathionine and maltotriose content increase in phase 3 13 . Alfalfa appears to cope with temperature increases mainly through decreases in various amino acid and glycerophospholipid contents, including tyrosine, succinic acid, methionine, glycerophosphate, glycerophosphocholine, lysophosphatidylcholine (16:0), and lysoPC (15:0) ( Table 2). These observations, combined with the results of previous studies by Xu 1,13 , highlight the importance of amino acids and glycerophospholipids in both cold acclimation and de-acclimation. Their www.nature.com/scientificreports/ contents increase during cold acclimation and decrease during de-acclimation. Soluble sugars are mainly involved in increasing the cold hardiness of alfalfa as part of the cold acclimation process.

Conclusions
Differences in the cold hardiness of two alfalfa cultivars gradually decreased in the continuous cold acclimation period, but increased during de-acclimation. These results indicate that the rapid decrease in the cold hardiness of plants in early spring is an important factor related to the inhibition of alfalfa regrowth. Results of the metabolomic analysis showed that ZD had better cold resistance than WL440 during both cold acclimation and de-acclimation, which was related to the change rate of amino acid and glycerophospholipids content in the two periods. In contrast, the primary function of soluble sugars is to improve the cold hardiness of alfalfa during the cold acclimation period. Data provided herein proves that soluble sugars, amino acids and lipids are important metabolites affecting the cold resistance of alfalfa, and it is further clarified that unsaturated fatty acids are conducive to improving the cold resistance. Furthermore, the study results provide a reference for breeding new cold-resistant alfalfa cultivars.

Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].